Experimental Implementation of Quantum Random Walk Algorithm 
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The quantum random walk is a possible approach to construct new quantum algorithms. Several 
groups have investigated the quantum random walk and experimental schemes were proposed. In 
this paper we present the experimental implementation of the quantum random walk algorithm on 
a nuclear magnetic resonance quantum computer. We observe that the quantum walk is in sharp 
contrast to its classical counterpart. In particular, the properties of the quantum walk strongly 
depends on the quantum entanglement. 

PACS numbers: 



I. INTRODUCTION 

Since the discovery of the first two quantum al- 
gorithms, Shor's factoring algorithm and Grover's 
database search algorithm 'J], research in the new born 
field of quantum computation exploded l^. Mainly mo- 
tivated by the idea that a computational device based 
on quantum mechanics might be (possibly exponentially) 
more powerful than a classical one great effort has 
been done to investigate new quantum algorithms and, 
more importantly, to experimentally construct a univer- 
sal quantum computer. However, finding quantum al- 
gorithms is a difficult task. It has also been proved ex- 
tremely difficult to experimentally construct a quantum 
computer, as well as to carry out quantum computations. 

Recently, several groups have investigated quantum 
analogues of random walk algorithms |E IE IE IE IE ^| j 
as a possible direction of research to adapt known clas- 
sical algorithms to the quantum case. Random walks on 
graphs play an essential role in various fields of natu- 
ral science, ranging from astronomy, solid-state physics, 
polymer chemistry, and biology to mathematics and com- 
puter science |Tlj. Current investigations show that 
quantum random walks have remar^bly__different fea- 
tures to the classical counterparts 5, 6, 7, 8, 9, 10]. The 
hope is that a quantum version of the random walk might 
lead to applications unavailable classically, and to con- 
struct new quantum algorithms. Indeed, the first quan- 
tum algorithms based on quantum walks with remark- 
able speedup have been reported ^^l- Further, experi- 
mental schemes have also been proposed to implement 
such quantum random walks by using an ion trap quan- 
tum computer Q and by using neutral atoms trapped 
in optical lattices Up to now, only three methods 

have been used to demonstrate quantum logic al g ates: 
trapped ions 0, cavity QED Tl^ and NMR [l||. Of 
these methods, NMR has been the most successful with 
realizations of quantum teleportation ^^Ij quantum er- 



ror correction ^^li quantum simulation quantum 
algorithms 0, quantum games and others [IJ. In 
this paper, we present the experimental implementation 
of the quantum random walk algorithm on a two-qubit 
NMR quantum computer, and we believe it is the first ex- 
perimental implementation of such quantum algorithms. 

We consider continuous-time random walks (CTRW) 
rather than discrete-time random walks ^ , on a circle 
with four nodes. We show that the evolution of this quan- 
tum walk is periodic and reversible, and yields an exactly 
uniform probability distribution at certain time. While 
the classical CTRW is irreversible and only approximates 
the uniform distribution at infinite-time limit. Further, 
we experimentally implement this quantum walk on a 
two-qubit quantum computer, using a unitary operator 
which has the right "effective" Hamiltonian, with good 
agreements between theory and experiment. It is inter- 
estingly found that the property of the quantum walk 
strongly depends on the entanglement between the two 
qubits. The uniform distribution could be obtained only 
when the two qubits are maximally entangled. 



II. QUANTUM CTRW ON A CIRCLE 

The concept of quantum CTRW is proposed in Ref . |E] ■ 
On a circle with four nodes, we denote the set of nodes 
by {0,1,2,3}. Since the structure of the circle is peri- 
odic, only two nodes, (fc -1-1) mod 4 and (/c — l)mod4, 
are connected to node k {k = 0,1,2,3). In the classical 
CTRW, let 7 denote the jumping rate, which is a fixed, 
time-independent constant that represents the probabil- 
ity moving from a given node to a connected one per unit 
time. The generator matrix j6j of this walk can therefore 
be written as 
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Consider a particle walking classically on this circle, let 
(t) denote the probability of being at node k at time 
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FIG. 1: The probabilities of being at the four nodes in the 
classical version of this CTRW. The solid (dotted) line cor- 
responds to P(p (t) (P2 (t)). The dashed line corresponds to 
Pf (t) and Pf (t), since Pf (t) = Pf (t). 



FIG. 2: The probabilities of being at the four nodes in the 
quantum version of this CTRW. The solid (dotted) line cor- 
responds to P^ (t) {P2 [t)). The dashed line corresponds to 
P^ (t) and P.^^ (t), since P^ (t) = P^ (t). 



t (the superscript C denotes ClassicaF ) , then we have |^ (g)) = |0), from equation Q, we have 



dt 



Pt it) 
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This equation conserves the probability in the sense that 

A natural way to the quantum version of this CTRW 
is to construct a Hilbert space spanned by the four basis 
{|0) , |1} , |2} , |3)}, respectively corresponding to the four 
nodes {0,1,2,3}. The state of a particle walking quan- 
tum mechanically on this circle is then denoted by IV' (i)), 
which is generally a superposition of the four basis, rather 
than classical mixing of the probabilities of being at the 
four nodes. The generator matrix in classical walks is 
now treated as the Hamiltonian of the quantum evolu- 
tion 0. Therefore the Schrodinger equation of the state 



;V (fc|ij|z)(/|^(0). (3) 



If measuring at time t, we can obtain a certain probabil- 
ity distribution on the circle. The probability of being at 
node k is (t) — \ {k\ if) {t))^ (the superscript Q denotes 
" Quantum" ) , and the conservation of probability is obvi- 
ously guaranteed by the normahzation (ip {t)\ ip (t)) = 1. 

Let the walking particle start from node 0, it is then 
easy to find the probability of being at node k at any 
time t in both the classical walks and the quantum walks 
(by solving equations Q and respectively). The de- 
tailed calculation yields that for the classical walks the 
probabilities are 
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= cos" 7t |0) - siT^^^it |2) 

+ ^e-2^^*sin27t(|l) + |3)). 

Therefore the probabilities of the quantum walks are 
{t) = cos^ 7i 

P!^ {t) = sin^jt 



(5) 



Pi^it) = P^it) = lsm'2jt 



(6) 



The probabilities in equations are plotted in FIG. 
n (for the classical walks) and FIG. [21 (for the quantum 
walks) as functions of time t. From FIG. ^ and FIG. 
13 we can see striking differences between quantum and 
classical random walks. FIG.[5]shows that the evolution 
of the quantum CTRW on a circle with four node are 
essentially periodic with a period T = t^/j- The par- 
ticle walking quantum mechanically on this circle will 
definitely go back to its original position, and the evo- 
lution is reversible. It is also interesting to see that at 
time t — 7r/27 the probability distribution converges to 
node 2. These phenomena are due to the quantum inter- 
ference effects, which allows probability amplitudes from 
different paths to cancel each other. 

To measure how uniform a distribution is, an immedi- 
ate way is to use the total variation distance between the 
given distribution and the uniform distribution. In our 
case, the classical and quantum total variation distance 
as functions of time t are 



In the quantum case, the initial state of the particle is 



(7) 
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FIG. 3: The quantum and classical probability distributions 
as functions of time t. The solid line corresponds to A*^ (t), 
and the dashed line to A'^ {t), both in theory. The crosses 
corresponds to the experimental results of the quantum case. 



FIG. 01 depicts the dependence of A"^ [t) and A'3 (t) on 
time t. From FIG. 13 we can see that the classical ver- 
sion of this walking process approaches the uniform dis- 
tribution exponentially as time lapses. In contrast, the 
quantum process exhibits an oscillating behavior. An in- 
triguing property of this quantum random walk is that 
A'^ (n7r/47) = if n is odd, which means that the prob- 
ability distribution is exactly uniform at time t = 7r/47 
and its odd multiples. While the classical walk can 
never reach the exactly uniform distribution, only ap- 
proximates it at infinite-time limit. 

III. EXPERIMENTAL IMPLEMENTATION 

For the quantum CTRW on a circle with four node, 
the Hilbert space is 4-demensional. So it is natural to 
implement the quantum walks on a two-qubit quantum 
computer. The direct correspondence is to map the basis 
{|0), 1 1), 1 2), 1 3)} of the quantum CTRW into the four 
computationalbasis {|0)«)|0), |0)(8)|1), |1)(8)|0), |1)®|1}}. 
This mapping is in fact to rephrase the number of nodes 
by the binary number system. Therefore the Hamiltonian 
in equation can be written as 

H ^2"fI®I -j{I(E)ax + (7x(^(Jx), (9) 

where / and are the identity operator and the Pauli 
operator of a single qubit. The evolution operator of the 
two-qubit system is 

U (t) = e-'"' 

= e"^*''* exp [ijt {ux ® (Jx)] exp [i-ft (/ ® Ox)] ■ 

(10) 

And the state of a particle performing this quantum 
CTRW is 

\'4, [t)) = e-2*T* cos^ -ft |00) - sin'^-it |10) 

+ le-2^7tsin27i(|01) + |ll)). (11) 
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FIG. 4: The correlation between the quantum total variation 
distance A*^ (t) and the entanglement S(t). The line corre- 
sponds to theoretical calculation. The experimental results 
in different sets are shown in different shapes. The trian- 
gles are for n £ {1,2,3}, crosses for n £ {4,5,6}, boxes for 
n G {7, 8, 9} and circles for n £ {10, 11, 12}. 



It is interesting to investigate the relations between the 
distribution of the implemented quantum CTRW and the 
entanglement of the two-qubit state \i){t)). The entan- 
glement of the two-qubit state \i' {t)) in equation 
can be directly calculated by the Von. Neumann entropy 

S {t) = — cos^ 7t log2 (cos^ 7t) — sin "^^t logj (sin^ 7i) . 

The correlation between the quantum total variation 
distance A'^ [t) and the entanglement S [t) is illustrated 
in FIG. 21 From FIG. 01 we can see that if there is no 
entanglement between the two qubits (5 = 0), A'^ is at 
its maximal A*^ = |, which corresponds converging at 
node (or node 2). While if the two qubits are maxi- 
mally entangled {S = 1), A'^ = 0, which happens to be 
the situation that the walking particle is uniformly dis- 
tributed on the four nodes. Therefore, we can say that 
the quantum random walk algorithm is enhanced by the 
quantum entanglement involved. 

The quantum CTRW is implemented using our two- 
qubit NMR quantum computer. This computer uses 
a 0.5 milliliter, 200 millimolar sample of Carbon-13 la- 
beled chloroform (Cambridge Isotopes) in dg acetone. In 
a magnetic field, the two spin states of and ^'^C nuclei 
in the molecular can be described as four nodes of two 
qubits, while radio frequency (RF) fields and spin-spin 
couple constant J are used to implement quantum net- 
work of CTRW. Experimentally, we performe twelve sep- 
arate sets of experiments with various selection of time t 
which is distinguished by 7t = (n = (0, 1, 2, • • • , 12}). 
In the following, we replace the jumping rate 7 with tt J 
(J — 215Hz). In each set, the full process of the quan- 
tum CTRW is executed. We describe this experimental 
process as follows. 

Firstly, prepare effective pure state \tp (0)): The initial 
state in NMR is thermally equilibrium state peq = ^ll + 
rather than a true pure state |V'(0)). However, it is 
possible to create an effective pure state, which behaves 
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in an equivalent manner. This is implemented as 

Rl (^/3) - - i?i (7r/4) - r - i?i (-7r/4) - G„ 

to be read from left to right, radio- frequency pulses are 
indicated by i?^^*"^ (angle), and are applied to the spins 
in the superscript, along the axis in the subscript, by 
the angle in the brackets. For example, R]. (tt/S) denotes 
7r/3 selective pulse that acts on the first qubit about x, 
and so forth. Gz is the pulsed field gradient along the z 
axis to annihilate transverse magnetizations, dashes are 
for readability only, and r represents a time interval of 
1/ (2J). Therefore, after the state preparation, we obatin 
effective pure state p (0) = ll + II + 21}; from equihb- 
rium state peq = 4/^ + . 

Secondly, perform quantum CTRW with different time 
t: As shown above in equation H10|l. quantum CTRW can 
be described as unitary operator U (t), this is performed 
with pulse sequence shown in the following (Note that 
the global phase e"^*''* of U (t) is safely ignored in our 
experiments, since p{t) = U (t) p (0) U (t)'' , this global 
phase is meaningless and has no effect on the result of 
experiment) 

Rl (6) - Rl' (^/2) - i - (^) „ ^ „ (_^/2) . 

Here i?^ (6) is equal to e"'^"^^ that act on the second 
spin, where 9 = nn/Q and t = n/ {6J) = nn/ (67) for 
n e {1,2, - ■ ■ ,12}, i?^^ (tt) denotes tt non-selective pulse 
that acts on both qubits about x. It is obviously that the 
final state p (t) of the quantum CTRW prior to measure 

is given by p{t) = U {t) p (0) U {t)\ 

Finally, readout the result p [t) and calculate quan- 
tum total variation distance A'^ (t); In NMR experiment, 
it is not practical to determine the final state directly, 
but an equivalent measurement can be made by so-called 
quantum state tomography to recover the density ma- 
trix p{t) = \ip {t)) {ip {t)\. However, as only the diag- 
onal elements of the final density operators are needed 
in our experiments, the readout procedure is simplified 
by applying gradient pulse before readout pulse to can- 
cel the off-diagonal elements. Here we shall note that, 
the gradient pulse can remove off-diagonal terms since 
we use heteronuclear systems in our experiment. Then 
quantum total variation distance A'^ it) is determined 



by the equation ls9 (t) = ^Yfk 



where 
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2 ^k=0 

(t) ~ {k\ p{t) \k) is certain probability on the node 
\k). Finally, A'3 (5") is determined with equation \i2\ . 

All experiments are conducted at room temperature 
and pressure on Bruker Avance DMX-500 spectrometer 
in Laboratory of Structure Biology, University of Science 



and Technology of China. FIG. El show the quantum to- 
tal variation distance A'^ (t) as a function of time t and 
FIG.^show the quantum total variation distance A'^ (S) 
as a function of entanglement S oi\il) {t)) shown in equa- 
tion (inj- From FIG. Hand FICH it is clearly seen the 
good agreement between theory and experiment. How- 
ever, there exist small errors which increase when time 
t increase, we think that the most errors are primarily 
due to decoherence, because the time used to implement 
quantum CTRW U (t) is increased from several to several 
tens milliseconds approximately, while the decoherence 
time T2 w 0.3 and 0.4 seconds for carbon and proton re- 
spectively. The other errors are due to inhomogeneity of 
magnetic field, imperfect pulses, and the variability over 
time of the measurement process. 



IV. CONCLUSION 

We present the experimental implementation of the 
quantum random walk algorithm on a two-qubit nuclear 
magnetic resonance quantum computer. For the quan- 
tum CTRW on a circle with four nodes, we observe that 
the quantum walk behaves greatly differently from its 
classical version. The quantum CTRW can yield an ex- 
actly uniform distribution, and is reversible and periodic, 
while the classical walk is essentially dissipative. Further, 
we find that the property of this quantum walk strongly 
depends on the quantum entanglement between the two 
qubits. The uniform distribution could be obtained only 
when the two qubits are maximally entangled. In this 
paper only the relatively simple case with two qubits are 
considered. However, our scheme could be extended to 
the case of a graph containing arbitrary N nodes, and 
the quantum random walk could be carried out by using 
log2 N qubits. 
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